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We analyze the phase structure of d = 4 A/" = 2 large N SYM theory with 
flavor on x by using geometric entropy as an order parameter. We introduce 
chemical potential conjugate to global U{1) symmetry and find the third order phase 
transition at finite density by using the geometric entropy as the order parameter. 
We also find that the geometric entropy at the finite density has an interesting 
behavior at low temperature and for large Nf. 
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1 Introduction 



Geometric entropy [T] was defined as the von-Neumann (information) entropy with Hamil- 
tonian associated with the space translation. The definition of this geometric entropy was 
shghtly different from that of the entanglement entropy [21 [3l H] which appeared in the 
condensed matter physics: geometric entropy was related to the entanglement entropy by 
double Wick rotation. However, geometric entropy becomes more convenient observable 
to investigate finite temperature systems including the strongly coupled Quark Gluon 
Plasma (QGP) and the weakly coupled deconfined phase of Yang-Mills theories. And 
particular interest is to investigate the phase structure of large QCD at a finite density 
by using geometric entropy as an order parameter. 

In the context of the gauge/gravity correspondence El E], geometric entropy was 
introduced as a new order parameter to analyze confinement /deconfinement transition of 
d = 4 large N SU{N) Yang- Mills theory on at finite temperature [1]. In the gauge 
theory side, free Yang-Mills theory [HI |9l [10] to mimic QCD has a natural scale, the radius 
R of S^, which is small enough to believe asymptotic freedom: RAqcd "C 1. And the 
confinement /deconfinement phase transition takes place if RAqcd is of order 1. 

In the dual gravity side, geometric entropy of the dual strongly coupled Yang-Mills 
theorjl^ obeys the area law formula. Namely, it is proportional to a minimal surface 
stretching on the AdS background. The first order confinement/deconfinement phase 
transition could also be seen holographically by using geometric entropy as an order 
parameter and the phase structure was consistent with the gauge theory results (see also 
the holographic calculation [I5]-[T7] of the entanglement entropy based on the AdS/CFT 
correspondence) . 

In this paper, we analyze the phase structure of free M = 2 large A^ Super Yang-Mills 
(SYM) theory with flavor on x by using geometric entropy as an order parame- 
ter. It is known that free M = 2 SYM with flavor can approximately be reduced to the 
Gross- Witten type matrix model [18], which can describe the third-order phase transi- 
tion [19]- [22] instead of the first order phase transition of the large A^ gauge theory. We 
also introduce chemical potential conjugate to global f/(l) symmetry and investigate the 
phase structure at a finite density. It is interesting to consider the application to the 
gauge/gravity correspondence [25] -[26]. However, M = 2 SYM theory is in general not 
conformal and so the application seems not to be straightforward except the conformal 
case A^/ = 2A^ or < A^. 

-•^ Interestingly, lattice theory results and the soft- wall AdS/QCD results [IT]-[T5 show that QCD is 
more like strong coupling. 
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The content of this paper is as follows. In section 2, we have shortly the review of 
geometric entropy. In section 3, we briefly review the partition function of free M = 2 
SYM theory with flavor on S*^ x S^. And we discuss the phase structure of this theory 
by obtaining free energy and Polyakov loop vev. In section 4, we capture the third order 
phase transition of A/" = 2 SYM theory with flavor by using geometric entropy. In section 
5, we introduce chemical potential conjugate to two global U{1) groups and generalize 
our model to the finite density system. 

2 Geometric entropy 




Figure 1: In the figure, we see the region near the orbifold fixed point 9 = tt/2. The 
deficit angle (a) is 5 = 27r(l — 1/ri). After a rescaling, there is the conical singularity (b) 
at ^ = 7r/2. 

Geometric entropy is an order parameter of the confinement /deconfinement phase 
transition for the finite temperature Yang-Mills theory on the compactified space. To 
describe the geometric entropy, we define the partition function on an orbifold S'^/Z„ 
to be ZYAiin): the ordinary partition function on x coincides with Zym{^)- The 
metric of in the unit radius i? = 1 is given by 

dQl = de^ + sin^ ed(f)'^ + cos2 ed-^^, (2.1) 

where < < 7r/2 and < 0, '?/^ < 27r. We define the Z„ orbifold action as follows: 

■0~?/^H . (2.2) 

n 

For n 7^ 1, we have the conical singularities at Q = n/2, with the deficit angle 5 = 
27r(l — 1/n) (see Figure. 1). The submanifold defined by these singular points is S^. 
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Geometric entropy is described by 

d 



G 



d{lln] 



log 



ZYM{.n) 



(2.3) 



n=l 



The geometric entropy is an analogue of the von-Neumann entropy as follows: 

ZYhiin) 



1)^ 



Trp", 5*0 = -Trplogp, 



(2.4) 



where p = e~ is the density matrix when we regard the coordinate (p as the Euclidean 
time and regard H as its Hamiltonian. Note that the geometric entropy is not equal to 
von-Neumann entropy using the Boltzmann distribution in a canonical ensemble. 



3 Free energy oi d = A gauge theories on an orbifold 

syz„ 

In this section, we compute the free energy of large U{N) gauge theories with flavor 
on the orbifold x S^fLn with the time direction compactified (t ~ t + /3) and with zero 
't Hooft coupling limit A = 0. We also calculate temperature dependence of the Polyakov 
loop vev in these theories. 

3.1 Partition function 

The matter contents of the d = 4 = 2 SYM are M = 2 vector-multiplet and M = 2 
hyper-multiplet: the hyper-multiplet contains 4 real scalars and 2 Weyl fermions. For 
free SYM theory on S*^, the matter field can be integrated out and the only remains are 
the temporal gauge field independent of positions. An exact expression for the partition 
function of the free M = 2 SYM on x 5*^ is described by the following unitary matrix 
model EOlET]: 



/[rf?7]exp ^^^^-^(^t;(x'")Trf/™Trf/™t + iAr/(x™)(Tr?7™ + Trt/'"t)j , (3.5) 

where / = {Nf/N)f' and Nf is the number of the fundamental-matters and Nf/N = 
const, and the sum over m denotes all the excitations of the single-particle partition 
function of the oscillators for adjoint v{x"^) and fundamental representations f'{x"^). 
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In the M = 2 SYM with flavor on the orbifold S'^ /'Ln-, assuming n is an odd integer, 
single-particle partition function of the matters for adjoint and fundamental representa- 
tion P, |24] are given by 

= vb^x"^) + {-iT+^VFix"^), (3.6) 
t^fiX = 2 — ^ o ^ + 2 ^ ^ , 3.7 

rrn/2+l 

vf{x) = 8- ^- -, (3.8) 

(1 — X) (1 — x") 

/'(x™) = /^(x-) + (-l)-+7].(x'"), (3.9) 

/^(a;)=4 ' /;(a;) = 8 ^2 ' (3-10) 

^ (l-x)'(l-x") ' (l-x)'(l-x") ^ ^ 

where x = e~'^ = e~^/^ at the temperature T and the index B, F means that we take the 
trace over the bosonic states and fermionic states, respectively. T is defined in units of 
the inverse radius l/i? to be dimensionless. 

We review these formulas on the orbifold S'^fZn following the appendix B of [9]. Since 
the states of the field theory on M x 5^ are related with the local operators on by a 
conformal transformation, the energy of the states is equal to the conformal dimension of 
the local operators. For this purpose, we embed the sphere in with the coordinates 
{zi,Z2)- Then, we operate the following Z„ orbifold: 

Zl ~ e'^zi. (3.11) 

For example, Z„ action operates on the scalar operator on as follows: 

0(2:1, Zl, Z2, Z2) ~ 0(e*^zi, e-*^zi, 22, 22), (3.12) 
9^0(2:1, Zl, Z2, Z2) ~ e*^9f0(e*^2;i, C^zi, Z2, Z2), (3.13) 

where we defined the partial derivative di = d/dzi {i = 1, 2). 

To obtain the single-particle partition function of the scalar fields, we sum the invariant 
operators under the Z„ action and the higher derivatives of the scalar operators with the 
traceless condition (the equation of motion). For n = 3, we should count the following 
operators: 

0, (920, ^20, <92<920, <92<920, <92<920, 5i5i9i0, dididi(f), ... (3.14) 

Note that the local operator ...didicj) is equivalent to ...92^20 because of the equation of 
motion. 
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The single-particle partition function is defined by 

z{x) = x^. (3.15) 

local operators 

When we use the book keeping parameter y and for di and di, respectively, the 
partition function of a real boson f l3.15p becomes 

x(l + x") 



(1 -X)2(l - X'^)' 

For the gauge fields on M^, the orbifold operates them as 



(3.16) 



A^-^{zi,zi,Z2,Z2) ^ e'~A^-,{e'~zi,e~'~zi,Z2,Z2), (3.17) 

- 27r ■ 27V 

A;,^{zi,zi,Z2,Z2) A^^{e'-zi,e-'—zi,Z2,Z2). (3.18) 



The single-particle partition function for is obtained by summing the operators 
obeying the identification didiA^ = —d2d2A^ and satisfying the condition, 

A'^ = 0, d^^.A^,) = 0, d^^,d,,A^,) = 0, (3.19) 

where () means symmetrization. These conditions are obtained by mapping the Gauss law 
constraint = on M x S''^ to x'^A^ = on and operating derivatives on x'^A^ = 0. 
For n = 3, we count the following invariant operators: 

diA,,, d2A,,, d2diA,,, d2d,A,,, ^s^aA,,, 92^2^,,, dfA,,, Of A,,,... (3.20) 

The single-particle partition function becomes 



z,{x) = 2x Efc=i kx'^il + 2 E,=i a;"'-^) 
_ 2x^(1 + 2x^^-1 -x") 
~ (1 -x)2(l -x'^) 

For Weyl fermions, the orbifold operates them including holonomy as 



(3.21) 



i7v{n-\~l) ■ 27r ■ 27v 

ip{zi, zi, Z2, Z2) e^^'^='^/;(e'~2i, e"*~^i, 22, ^2), (3.22) 

where Weyl fermions are of 2-component and operated by (T3. Note that the fermions are 
rotated by 1/2 of bosonic variables. Since every fermionic field is affected by the orbifold 
action, the invariant operators for n = 3 are of the following form: 

diipi, dii)2, 8281^)1, d2diipi,.... (3.23) 
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The single-particle partition function of a Weyl fermion becomes 



Ax 



x) 



X" 



(3.24) 



where 4 in the numerator of fl3.24p comes from contributions of complex fermions tjji and 

^2- 



3.2 Free energy 

We derive the free energy of A/" = 2 SYM by using the approximation v{x^) = f{x^) = 
(m > 1). Namely, the first winding state in the time direction is only excited. This 
approximation is valid for not sufficiently high temperature region. By introducing the 
Lagrange multiplier A, A [HI [23], (13. 5p is rewritten as 

Ar2 



Z{vJ) 



inv 
Ar2 



4:71V 



f[dU]dXd\ exp 



- N - 

— {X-m-f) + -{XTrU + XTTW) 



- — {g^-2gfcos9 + p) 



Ng, 



■ J[dU] exp l^-yiTiU + Trf/t) 
2v ■'^ 



iV2 



— C gdge 



-N^fiF{v,f,9) 



where 

N^I3F{vJ,g) 
and we defined 



j [dU] exp ]^Ng {TiU + Trf/^) . 



(3.25) 
(3.26) 
(3.27) 
(3.28) 

(3.29) 
(3.30) 



The asymptotic expansion of K{g) in the large limit is computed in [18] as follows: 



g--\ogg-- + 0{l/N^) {g > I) 
We also expand the first term in fl3.29p in powers of A^^ as follows: 

'N^gf 



(3.31) 



log/o 



2v 



2v 2 V V 



m^gf 



0{N- 



(3.32) 
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I Nf;N=0.5 Nf;N=0.1 Nf;N=2| 



Figure 2: Temperature dependence of the free energy for the J\f = 2 SYM with flavor is 
plotted. 

We analyze the phase structure by using the saddle point approximation and by looking 
at the leading order terms in (13.291) in powers of N"^. 

When f < 1 and f < fo = 1 — v , F{v, f, g) has a minimum at 

^7o = T<l' PFi^J^9o) = - ,J' y (3.33) 
/o 4(1 - V) 

This behavior of free energy implies that the fermions are not the fundamental degrees 
of freedom but the mesons are fundamental degree of freedom N'^F oc Nj at the low 
temperature. 

When V < 1 and / > /o or f > 1, F{v, /, g) has a minimum at 



= » + | + W„./.,.) = -|-/f + i + ilog9„ + £.(3.34) 

We plot T-dependence of the free energy for various A'^^ in Fig. |2J 

Note that the critical line of the above phases is described by the following formula: 

v{x) + fix) = 1, go = 1. (3.35) 



As was shown in the literature [181 120], there is the Gross- Witten type third order 
phase transition at the temperature T = Tc determined by the formula (I3.35p . 
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T 

Nf/N=0.5 Nf;N=0.1 Nf/N^ 



Figure 3: Temperature dependence of Polyakov loop vev is plotted. T is in units of 1/R. 
The critical temperature Tc at which P = 1/2 is equal to 0.46,0.41,0.34 for Nf/N = 
0.1, 0.5, 2, respectively. 

3.3 Polyakov loop vev 

It is interesting to obtain the Polyakov loop expectation value. According to [19], the 
density of eigenvalues becomes in our notation, 

{g a (I .2 a\\ • 2 a / i 

— cos — sm — 9 > 1, sm 77 < -, 
— (l+5(cosa) g<l, |a| < vr. 

Polyakov loop vev becomes the first moment of p as follows: 

pcos{a)da. (3.37) 

-TT 

And it can be obtained from f l3.36p as 

p, = { pWKlf^ (3.38) 

^ < 1- 



2(1 -t; 



We plotted temperature dependence of Polyakov loop vev for Nf/N = 2, 0.5, 0.1 in Fig. [3l 
It is interesting to compare Fig. 12] with the Figure in [12] . Our results of Polyakov loop 
vev have a similar behavior with the lattice results and the soft-wall AdS/QCD results 
written in [T2]. 
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4 Geometric entropy of = 2 SYM 




I Nf/N=0.5 Nf/N=0.1 Nf/N=2| 



Figure 4: Temperature dependence of geometric entropy is plotted. T is in units of l/i?. 




-30- 

I Nf/N=0.5 Nf/N=0.1 Nf/N=2| 

Figure 5: Temperature dependence of dS^ /dT is plotted. T is in units of 1/-R. The phase 
transition happens at = 0.46, 0.41, 0.34, for Nf/N = 0.1, 0.5, 2, respectively. 



By using the free energy f l3.33p and f l3.34p . we can describe geometric entropy of A/" = 2 
SYM as follows: 



Sr 



d 



d{l/n 

_d_ 

dn 



\ogZ{n)--\ogZ{r 
n 



(/3F)(n)-i(/3F(l))). 



n=l 



(4.39) 
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where Z{n) = Z{v, f) and F{n) = F{v, f, go). Computing fl4.39p . we plot Sg and dSc/dT 
near the critical temperature in Figure. H] and Figure. |5l respectively. We find that 
since geometric entropy is a kind of von-Neumann entropy, geometric entropy of A/" = 2 
SYM can capture the third order phase transition of the Gross- Witten matrix model and 
the critical temperature Tc decreases if Nf/N increases. Tc is consistent with the value 
obtained in fl3.35p . 

It will be interesting to consider the case at non-zero 't Hooft coupling. According 
to [21], the order of the phase transition will become the first order at non-zero 't Hooft 
couping. So, it is expected that in this case, the geometric entropy can capture the first 
order confinement /deconfinement transition. 



4.1 Analysis for the high temperature Umit 

In the high temperature limit a; — ?■ 1 (/3 — j- 0), the single-particle partition functions on 
S^/Zn dMU) and (13:2^1) are given by 



/ m\ 4 - 6n - 1 

^ = + + ^ ^ ' 

4 2 -I- r?2 

zr(x"') = ——--—J- + 0(f3). (4.40) 

where Zc{x) is the single-particle partition function of a complex boson. Using (14. 40 p . the 
free energy and the geometric entropy of the M = 2 SYM are obtained as follows (see 

my 

^ = -^^^iV^r^(l + ^)v^53, (4.41) 



6/3 V ^ 

where 1^53 = 27r^. The above result implies that the geometric entropy becomes zero 
similar to the case of A/" = 4 SYM if = Nf and the geometric entropy changes its sign 
depending on A^ and A^/@ 



^The coefRcient in (|4.42p is similar to the /3 function of A/" = 2 SYM theory. However, in the QCD 
case with no supersymmetry, the geometric entropy has minus value for any Nf since the only positive 
contribution to the geometric entropy is obtained from that of scalars and periodic fermions (see [1]). 
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According to [T], So should be subtracted by the geometric entropy Sq, where Sq is 
geometric entropy with the fermions obeying the periodic boundary condition in the time 
direction: tr(— l)^e~'^^ instead of tre"'^^. The partition function with periodic fermions 
is easily obtained by removing (—1)'^ in (13.61) and (13. 9p . Thus, Sq and ASq = Sq — Sq 
is obtained as follows: 

SS-q^. (4.43, 



TT 



'■'N'' (N, 



The result (I4.44p is plausible since the coefficient of AS'g is proportional to that of the 
entropy —dF/dT. 

We want to comment on the relation to the gauge/gravity correspondence for Nf/N ^ 
1. ASg seems to be comparable with the dual gravity result in the high temperature limit. 
However, the analysis in the dual gravity side seems to be difficult since in the probe D7- 
brane analysis [25| [26] , geometric entropy has no flavor contributions. It implies that 
higher order correction on the D7-brane worldvolume and the background charges of the 
D7-brane should be included. 



5 Introducing chemical potential 

In this section, we analyze the phase structure of c? = 4 A/" = 2 SYM with flavor at flnite 
chemical potential. The similar analysis in d = 4 M' = 4 SYM was done in |271 EHl [29] . 
We start with reviewing global symmetry oi d = 4 N' = 2 SYM with flavor. The form 
of the Lagrangian for Nf hypermultiplets {Qa{q,i^q),Qa{Q,'ipq)) (labeled by an index a), 
interacting with a A/" = 2 vector multiplet is given by 

C = J (Qle-2^g, + Oae'^gl) + J d''9iQa<^Qa) + h.c. (5.45) 

where V is the A/" = 1 superfleld and $ is chiral superfield in the adjoint representation. In 
M = 2 SYM theory, there are U(l)j and U{1)r subgroup^ in SU{2)r i?-symmetry [30] . 
In this paper, we don't consider chemical potential conjugate to U{1)r since M = 2 SYM 
with flavor on 5^ x S'^ will be not free from the anomaly for large Nj. The transformations 
of A/" = 1 superflelds under U{l)j are given by 

U{l)j ■ $ ^ <I)(e-*°^), V V{e-''^e), (5.46) 

Q e^°g(e-^"^), Q e*"Q(e-^"^), (5.47) 

■^Remind that U{l)ii operates on A/" = 2 fields as a chiral symmetry that is quantum- mechanically 
broken to the discrete Z^m-^Nj due to anomalies except the case Nf = 2N . 
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And the transformations of all component fields under U{l)j are given by 



(A,g)^e-(A,g), 



(5.48) 



where A) are fermions in the M = 2 vector multiplet. 

In addition, there is an U{\)f subgroup in U{Nf) flavor symmetry as follows: 



U{1)f symmetry is baryonic symmetry of A/" = 2 SYM with flavor. 

In the next section and the appendix A, we introduce chemical potential conjugate to 
these f/(l) symmetry and analyze the phase structure oi M = 2 SYM at finite density. 



Figure 6: Temperature dependence of geometric entropy is plotted at finite chemical 
potential /i = 1/2. In the figure, T and n are measured in units of 1/R. Our plot shows 
that Sg at low temperature and for large Nj takes the positive value. 

For U{l)j case, we make the approximation of ignoring the higher order of winding 
states m > 1. This approximation is valid for not high temperature and not at large 
densit}!^. Then, the phase structure of the system at finite chemical potential can be 
analyzed by evaluating free energy, geometric entropy, and number density. 

^Note that there are also stabihty bounds E > ^Q, where E is the energy of particle states and Q is 
an U{1) charge of these states: in (|5.54p . the bound is < 1. 



(g„gl)^e^"(g„Ql). 



(5.49) 



5.1 U{l)j case 




0.0- 
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-20- 

I Nf;N=0.1 Nf;N=0.5 Nf;N=2] 



Figure 7: Temperature dependence of dSc/dT is plotted at finite chemical potential 
/i = 1/2. T and fi are in units of 1/R. The phase transition happens at Tc = 0.41, 0.34, 
and 0.22 for Nf/N = 0.1, 0.5, and 2, respectively. 



0.2 "0:4- 



0.6 Tr.S^^ 1.0^^[.2 
' ■ - -mu \ 



\ 



\ 



I Nf/N=0.1 Nf;N=0.5 Nf;N=2| 



Figure 8: /i dependence of F at T = 0.3 is plotted, fi is measured in units of 1/R. 

For single-particle partition functions including chemical potential fi conjugate to the 
U{l)j charge, we should replace f l3.6p - fl3.10p by 

v{x"') = VBix"") + i-l)"'+hF{x"'), (5.50) 

..(X) ^ 2 + + 2 , (5.51) 

^ ^ (l-x)'(l-x") (l-x)'(l-x") ^ ^ 
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Figure 9: /i dependence of Sq sA, T = 0.3 is plotted, /i is in units of 1/R. Our plot shows 
that Sg at low temperature and for large Nf takes the positive value as /i becomes large. 
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I Nf/N=0.1 Nf/N=0.5 Nf/N=2| 



Figure 10: fi dependence of (PSc/d^"^ at T = 0.3 is plotted, fi is in units of 1/R. The 
phase transition happens at the critical values of chemical potential fic = 0.89, 0.62, and 
0.23 for Nf/N = 0.1, 0.5, and 2, respectively. 

n/2+l+n n./2+l-M 

vf{x)=4: 5 +4 5 , (5.52) 

^ ^ (l-x)'(l-x") (l-x)'(l-x-) ^ ^ 

f (x™) = /^(a;™) + (-l)"+7;.(x"), (5.53) 

/^(a;)=4^^ \ ^ , /^(x) = 8 . (5.54) 

' (1 - a:)' (1 - X") ^ (l-x)'(l-x") ^ ^ 

Here, we introduce the number density b depending on fi as an order parameter of the 
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Figure 11: /i dependence of the number density b is plotted at finite density and T = 0.3. 
fj, is in units of 1/i?. b increases monotonously as fj, increases as expected. 

200- 




- 100-1 

I Nf;N=0.1 Nf;N=0.5 Nf;N=2| 



Figure 12: /i dependence of d'^b/dfi'^ is plotted at finite density and T = 0.3. fi is in units 
of 1/R. 

third order phase transition. We define b as 

b=-^. (5.55) 
djj 

The free energy can be evaluated by calculating (13.331) and (I3.34p in the presence of 
chemical potential and has behavior similar to /i = for 7^ 0. We plotted temperature 
dependence of S'g at /i = 1/2 and dSc/dT at fi = 1/2 in Fig. [HI and Fig. O respectively. 
We observe the third order phase transition in the presence of finite chemical potential 
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/i=l/2. 

The geometric entropy for the large flavor Nj/N > 1 has interesting behavior at low 
temperature. For the large flavor Nf/N = 2, the geometric entropy takes the larger 
positive value at low temperature than Nf/N = 0.1, 0.5 cases. The similar behavior is 
also observed evaluating Sq and dSc/dfi as the functions of /i in Fig. [9] and Fig.[TOl These 
behavior of the geometric entropy at the low temperature seems to come from the relation 
to the scale anomaly that is proportional to the central charges. Since the definition of 
the geometric entropy includes changes of the area Area(5''^/Zi„) by varying the orbifold 
action, the geometric entropy at the low temperature will be related with a scale anomaly 
of AT = 2 SYM theory on M x 

Moreover, /i dependence of b and (Pb/dfi^ at T = 0.3 is plotted in Fig. [11] and Fig. [T2l 
respectively. Our results show that both b and Sq are order parameters of the third order 
phase transition. 



6 Discussion 

In this paper, we investigated the phase structure of c? = 4 A/" = 2 large SYM theory 
with flavor at finite temperature by using geometric entropy as an order parameter. Ana- 
lyzing geometric entropy, namely von-Neumann entropy associated with the space trans- 
lation, we observed the third order phase transition as the confinement /deconfinement 
transition instead of the first order phase transition. We also enlarged our analysis to 
the finite density system and observed the third order phase transition in respect to T 
or /i. The geometric entropy at finite chemical potential for large Nf has the interesting 
behavior at low temperature implying that the scale anomaly of A/" = 2 SYM on M x S*^ 
may contribute. 

In the appendix A, we also analyzed the system in the presence of the baryon chemical 
potential conjugate to U{1)f in the low temperature limit since the similar analysis at 
finite temperature is more complicated. We found the confinement phase similar to the 
quark matter phase in the large A^ QCD at the finite chemical potential n > 1. For future 
works, we should include the interaction terms of bosons to make the baryon number 
density with positive chemical potential /i > well defined. 

It will be interesting to include mass terms for the flavor field as had been done in [22] 
for large A^ and A^ = 3 QCD theory since the (T, /i) phase space of two theories has 
different behavior. Moreover, it is important to compute the chiral condensate of A/" = 2 

^Thc period of is large at low temperature: ~ M. 
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SYM theory with massive flavor fields and to argue the chiral symmetry breaking of this 
theory. It will also be interesting to enlarge our analysis for different orbifold systems [21] 
since geometric entropy depends on the orbifold singularity. 
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A U{1)f case and low temperature limit 

U{1)f symmetry is similar to the baryonic symmetry in QCD theories. For the U{1)f 
case, it seems to be interesting to investigate finite density and the low temperature phase 
of SU{N) gauge theories such as quark matter phase in the large N QCD [32], while 
the analysis at finite temperature by using geometric entropy is complicated. In this 
subsection, the (T, /i) phase space near /i = 1 and at low temperature is investigated by 
using saddle point approximation. Introducing chemical potential conjugate to U {1)f, the 
chiral field Qa in the M = 2 hypermultiplet survives at low temperature (-R//3 ^ 1) [22] . 
In addition, fermionic modes and higher order of scalar modes can be ignored in the small 
sphere limit. In the low temperature limit, the partition function becomes 

z = jY[,m^M-sm). (A.56) 

S{ed = -^E.^logsin^ (^^) +iVE.^(^.)> (A.57) 
V{e) = ai log(l - ^e*^) + iAfe, (A.58) 

where ^ = e'^'-^"^-' and we introduced the Lagrange multiplier enforcing the ^ ■ 9i = 
constraint and focused on the lowest energy level of the bosonic states, ai = 2Nf/N% 

Introducing the parameter s satisfying —it < s < it and s G M and defining Zi = 
exp{i6i), the eigenvalue distribution in the large N limit is defined by 

t^=p{z). (A59) 
az 



^The generalization to the multi-level model is straightforward. However, it will be more complicated 
to find the boundary in the {fi, T) plane. 
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The saddle point equation thus becomes 

zV'{z) = pf^p{z'f-±^, zV\z)=M-^. {km) 
J Zm z — z' i- — c,z 

p{z) has poles at 2; = and 1/^. There are two cases where the pole at 1/^ is outside the 
contour of integral or not. Note that these two cases are related with each other by the 
transformation ^ — 1/^. So, it is enough to consider the small ^ case where 1/^ is outside 
the contour. In the small C, limit, the pole at 1/.^ is outside the contour of integral. Then, 
the following p{z) satisfies flA.60P : 

p(.) = i + ^. (A^61) 

Polyakov loop vev is given by 

/dz 
— p(^)^ = 0. (A.62) 

Vanishing Polyakov vev implies that the confinement happens. The phase transitions to 
the deconfinement phase happen at the place where p{z) vanishes precisely at the point 
z{±tt). The line of the phase transitions in the (/i,T) plane is given by 

fi=l + T{{ai - 1) log(-l + cri) - logai). (A.63) 

To make (1A.63I) well defined, the condition Nf/N > 1/2 should be satisfied. And the 
confinement region at finite chemical potential exists if Nf/N > 1/2 implying that the 
confinement phase of A/" = 2 SQCD is derived at least by evaluating the fermionic states. 
The next confinement region is separated by the line 

p = l- T{{ai - 1) log(-l + ai) - ai logai). (A.64) 

The above line is obtained by replacing ^ with 1/^ in (1A.58I) . Note that the line ( 1A.63I) 
and flA.64p reach /i = 1 as Nf/N — )■ 1/2. And the deconfinement phase vanishes at the low 
temperature. This behavior is similar to that of the large QCD in a small analyzed 
in [22]. For these cases, it can be shown that in the small limit and as Nf/N — > 0, the 
corresponding matrix model can realize the phase structure of the large QCD in the 
low temperature limit This phase structure contains the quark matter phase where the 
pressure behaves as 0{N) and the Polyakov loop vev vanishes. 

^According to [35], furthermore, the baryon degeneracy at large N has the function similar to the 
coefHcient of T in (|X63l) . 
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